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MEASUREMENT OF MAGNETIC FIELDS USING A STRING FIXED AT BOTH ENDS 

5 The present invention relates to the measurement of 

magnetic fields, particularly to the measurement of 
components of the magnetic gradient tensor, the 
measurement of which is known as magnetic gradiometry. 
Accurate absolute measurements of the spatial 

10 derivatives of magnetic induction, i.e. the magnetic 

gradient tensor = dB ± /dj (ij = x,y,z), taken at some 
local coordinate frame OXYZ are very important to 
progress in the fields of geological prospecting, sea 
and underwater navigation and exploration, terrestrial 

15 and marine archaeology and medicine ( magnet ocardiography 
and magnetoencephalography) . One reason for this 
importance in relation to geology and archaeology is 
that the Earth's magnetic field is approximately uniform 
in most regions of interest, so that magnetic gradient 

20 measurements are unaffected by the Earth's magnetic 
field. 

There do not currently exist any devices which are 
capable of directly measuring magnetic gradients. One 
known indirect method uses two or more spaced 

25 magnetometers to take readings of the magnetic 

induction. The measurements from the magnetometers are 
differenced and divided by the spacing at which the 
readings were taken, to give an approximation of the 
magnetic gradient. The magnetometers are usually spaced 

3 0 by around 5 0cm, which leads to a relatively large 

apparatus. In addition, two magnetometers are required 
for each direction along which the gradient is to be 
measured, which increases the cost of the apparatus. 
Magnetic flux gradiometers exist which utilise 

35 SQUIDs (Superconducting Quantum Interference Devices) , 
but these require cryogenic temperatures and finely 
balanced superconducting pick-up coils, which limit 
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their practical usefulness and accuracy. 

The present invention utilises a novel sensitive 
element in the form of a current -carrying string and 
suitable sensors to detect magnetic fields, in 
5 particular magnetic gradients. 

Thus, viewed from a first aspect the invention 
provides apparatus for the measurement of magnetic 
fields comprising: 

a string fixed at both ends, held under tension and 
10 arranged to carry a current; 

and sensing means arranged to detect the transverse 
displacement of said string from an unperturbated 
position, due to a magnetic field acting on said string. 

By "string" no particular limitation as to material 
15 or construction is intended. Any elongate element is 
included which is capable of carrying a current, of 
being transversely deflected by a magnetic field and of 
providing a restoring force . 

Viewed from a further aspect, the invention 
2 0 provides a method of measuring magnetic fields 
comprising : 

providing a string fixed at both ends, held under 
tension and carrying a current; and 

detecting the transverse displacement of said 

25 string from an unperturbated position, due to a magnetic 
field acting on said string. 

An unperturbated flexible string with fixed ends 
forms an absolute line in space going through the points 
where the ends of the string are fixed. This line can 

30 be arbitrarily selected as one of the axes of a local 
cartesian coordinate frame, say Z, with the other two 
axes, X and Y, chosen to lie in the plane transverse to 
the direction of the string. Any string deflection from 
this line is caused by absolute values of the transverse 

35 components of the force per unit length which acts on 
each unit element of the string. The force per unit 
length is due to the ambient gravitational field and, 



10 
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when the string carries a current, to the ambient 
magnetic field. 

The displacement y(z, t) of a string of length 1 
from its unpertubated position, for example, in the y- 
direction of the above local coordinate frame as a 
function of the z-position of a unit element and time 
can be described by the following differential equation 

d 2 / u 3 , , , v , Al d 2 . 

ri -y(z, t) + h — y{z, t) - YA— -y(2, t) 

dt 2 ot 1 dz 2 

= - ng y (0, t) - nr yz (0, t) z 

+ I(t)B x (0,t) + I(t)B xz (0, t) z (1) 

+ noise 

15 with boundary conditions corresponding to the fixed ends 
of the string, i.e. y(0,t) = yd, t) = 0. In this 
equation r\ denotes the string's mass per unit length, h 
is the friction coefficient per unit length, the 
parameters Y, A and Al/2 are the string's Young modulus, 

20 the area of its cross-section and the string's strain 
respectively. The quantities g y (0,t) and r yz (0,t) are 
the absolute values of the y-component of the total 
acceleration and the corresponding gravity gradient 
tensor component along the string, both taken at the 

25 'centre of the local coordinate frame chosen. 

The quantity J (t) is the current flowing through 
the string. It is well known that a conductor carrying 
a current I(t) in a non-uniform magnetic field of flux 
density B(x,y, z) is subject to force F = I(t){n x 

30 B(x,y, z) } , where n is the unit vector in the direction 
of current flow, in this case the Z direction. The 
quantities B x (0,t) and B xz (0,t) therefore represent the 
absolute values of the x-component of the magnetic 
induction and the corresponding magnetic gradient tensor 

35 component along the string, both taken at the centre of 
the local coordinate frame chosen. 

In this example, the y-direction has been chosen as 
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10 



15 



an arbitrary example to simplify the explanation of the 
invention. However, the foregoing and following 
analysis is equally applicable to any direction 
transverse to the string or any number of directions. 

Applying Fourier analysis to the complex shape of 
the string caused by its interaction with the magnetic 
and gravitational fields, the function y(z,t), can be 
described, in the range z=0 to z=2, by an infinite sum 
of sinusoidal functions of period 21, with appropriate 
coefficients c y (n, t) . Thus, a solution of Equation(l), 
which satisfies the boundary conditions shown above, can. 
be represented by the following infinite sum wherein 
each term in n corresponds to one of the string's 
natural vibrational modes 



y(z,t) = ^ c y (n, t) sin^ iyzj 



(2) 



20 



By substituting Equation (2) into Equation (1) and 
by multiplying both sides by sin(nn'z/l) , and then by 
integrating both sides over z from 0 to 1, one can 
obtain the differential equation for c y (n,t) 



25 



d 2 / ^ 2 d . , x 2 , 
-c In, t) + - — c (n, t) + co„c (n, t) = 

dt T dt y y 



— [(-D n -i] 



nn 

nn 



g y (0, t) + -I(t)B x (0, t) 



r v2 (0,t) + -i(t)B(0, t) 

yz ^ xz 

+ noise 



(3) 



30 



where the quantities 



nn 



Y AI 
P 1 



(4) 



35 



represent the string's natural frequencies, and t and p 
are the relaxation time and the volume mass density of 



WO 00/68702 PCT/EP00/04270 

- 5 - 

the string respectively. Thus, when n=l, co n represents 
the fundamental resonant frequency of the string. 

The displacement of the string from its 
unperturbated position may be measured in the y (and or 
5 x) direction by one or more suitable sensors positioned 
along the length of the string. The deviation of the 
string from a straight line may be approximated in terms 
of a finite Fourier sum of the string's natural modes. 
For example, the deviation of the string from its 

10 unperturbated position may be recorded by an optical 
sensor or the like and analysed. According to one 
possibility, an image of the string in its deviated 
position may be recorded by a camera.. 

Alternatively, electrically inductive sensors may 

15 be used to detect the displacement of parts of the 

string. From this analysis, values for c y (n,t 0 ) may be 
calculated for a given time t 0f and in turn the values of 
any of g y (0,t 0 ), r yz (0,t 0 ), B x (0,t 0 ) and/or B xz (0, t 0 ) may 
be calculated. In this way, magnetic induction and 

20 magnetic gradient may be measured using the string. 

However, when n takes an even value, i.e. for those 
terms of the infinite sum in Equation (2) corresponding 
to vibrational modes of the string having a node at 
z=l/2 ( ant i- symmetric modes), the term involving g y {0,t) 

25 and B x (0,t) is equal to zero. When n is even therefore, 
c y is dependent only on r yz f0,t; and B xz (0,t) (and thermal 
noise) . 

Thus, by measuring only the amplitudes of 
antisymmetric modes in the Fourier sum, values for 

30 r yz (0,t 0 ) and B xz (0,t 0 j may be calculated. 

The sensing means may therefore comprise at least 
two sensors located symmetrically about the midpoint of 
the string. In this way, the sensors may be arranged to 
detect the amplitude of the ant i- symmetric natural modes 

35 of the string in its deflected position, for example by 
differencing the output from each sensor. In a 
preferred arrangement, the sensors are located at 
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positions corresponding to the antinodes of at least one 
anti- symmetric natural mode of- the string, for example 
n=2 , where the displacement of the string will be at a 
maximum . 

5 The current through the string may be an 

alternating current. In this case, according to 
Equation (1) the string will undergo forced oscillations 
generated by the force due to the alternating current in 
the magnetic field. The terms in B x and B X2 in Equation 

10 (1) will vary over time with the current I (t) , whereas 

the terms in g y and F yz will remain constant for a quasi- 
static gravitational field. This variation may be used 
to distinguish displacement of the string due to the 
magnetic field from that due to the gravitational field, 

15 and thereby quantify the magnetic field. 

Advantageously, the frequency of the alternating 
current may be selected to cause the string to resonate 
in one or more of its natural modes. In this case, the 
amplitude of the oscillation may be measured and will be 

20 equal to c y (n,t) if the string is resonating in the nth 

mode. As explained above, from the values of c y (n,t) the 
values of any of g y (0,t 0 ), r yz (0 f t 0 ), B x (0,t 0 ) and/or B X2 f0, 
t 0 ) may be calculated. 

In a particularly advantageous arrangement, the 

2 5 frequency of the alternating current may be 

substantially equal to the frequency of one of the anti- 
symmetric natural modes of the string, i.e. the 
frequency of the alternating current may be an even 
multiple (n) of the fundamental frequency of the string. 
30 For example, if the current through the string I(t) 

is selected to satisfy I(t)=i sin((x> n t), where co n is the 
frequency of the nth natural mode of the string (where n 
is even) , the string will undergo forced oscillations 
which will lead to resonance of the string. At 

3 5 resonance, when damping is completely counteracted by 

the pumping force of the alternating current, the 
amplitude of the oscillations will be given by the 
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following expression 

c(n,t) = — — —B cos (0) t) /,-x 
y n co n n n xz n 

5 Thus, the string resonates in the nth natural mode 

and according to Equation (5) the amplitude of the 
resonant oscillations is directly proportional to the 
magnetic gradient B xz , and can therefore be used to 
measure this quantity. For example, sensors located at 

10 the antinodes of this mode, i.e. symmetrically about the 
midpoint of the string, may be used to measure the 
amplitude directly. 

In a preferred embodiment, the frequency of the 
alternating current is twice the fundamental frequency 

15 of the string, i.e. n=2 . 

The string's deflection from its unperturbated 
position may be detected by any suitable displacement 
sensing device. When the string carries an alternating 
current, the sensing means may detect the displacement 

20 of the string by electromagnetic induction. The 

amplitude of the current induced in a conductor adjacent 
the string will be a direct function of the distance of 
the string from that conductor. Thus, in a preferred 
embodiment of the invention one or more fixed pick-up 

25 coils are arranged along the length of the string to act 
as displacement sensing means, the current induced in 
each coil being directly related to the string's 
displacement from its unperturbated position. 

In a particularly advantageous embodiment, 

30 displacement sensors, for example pick-up coils, are 
arranged adjacent the string in two non-parallel, 
preferably orthogonal, planes, so as to be capable of 
measuring the string's displacement in two transverse 
directions simultaneously. 

35 Advantageously, the alternating current may be used 

to demodulate the signals from the displacement sensors. 
For example, when the frequency of the alternating 
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current is an even multiple of the fundamental frequency 
of the string, the frequency of the signal generated in 
an electromagnetic pick-up element, for example a coil, 
adjacent the resonating string is a multiple of the 
5 signal due to the alternating current through the string 
in its unperturbated position, say sin (cd n t), with the 
signal due to the physical oscillation of the string, 
say -cos(u> n t), which is the integral of the alternating 
current function. This multiplication results in a 

10 signal from the sensors of frequency 2co n . A demodulating 
reference signal of frequency 2co n can be generated from 
the alternating current to synchronously detect the 
sensor signals. Because both signals are generated from 
the same source, it is ensured that the frequency 

15 relationship is constant. 

The voltage excited in the pick-up coils due to 
electromagnetic induction, is given by the following 
expression 

V e = GxB xz xcos(2a) n t) (6) 

20 

where G is the gain in the primary inductive transducer, 
for example a pair of pick-up coils. The value of G is 
defined by the geometrical configuration of the pick-up 
coils and is proportional to the current through the 

25 string. 

In a preferred embodiment, the string and 
preferably also the sensing means is cooled to reduce 
the effect of thermal noise on the measurements. For 
example the apparatus may be immersed in a bath of 

30 liquid nitrogen or other suitable coolant. 

The apparatus may comprise more than one string, 
for example three, four, five or more strings. The 
strings may be arranged at arbitrary angles to each 
other in order to measure simultaneously all the 

35 independent elements of the magnetic gradient tensor. 

Some embodiments of the invention will now be 
described by way of example only and with reference to 
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the accompanying drawings, in which: 

Figure 1 is a general schematic representation of a 
preferred embodiment of the invention; 

Figure 2 represents an experiment carried out using 
5 the embodiment of Figure 1 # and 

Figure 3 shows experimental results achieved with a 
second embodiment of the invention. 

A single channel prototype of a device according to 
the invention is shown in Figure 1. The device 
10 comprises a flexible string 1. The string has, in this 
embodiment, a length 1 = 26 cm, is 0.25 mm in diameter 
and is fixed at its ends in a ceramic housing. The 
string is preferably formed of a wire of low mass-per- 
unit length and relatively high tensile strength. In 
15 the example shown, the string is formed of composite 
zinc clad silver wire with a mass-per-unit length of 
0.44 x 10" 3 kgm" 1 and a tensile strength of 13 0, 000 psi 
(90 0 MNm" 2 ) . 

The conductive string is oriented vertically and 

20 two pick-up coils LI and L2 are connected in series and 
directly coupled t.o a low noise amplifier A. The pick- 
up coils LI, L2 are connected in antiphase such that 
identical signals which are induced in each coil cancel 
out, whereas opposite signals are added together. This 

25 configuration (when demodulation is carried out at a 

frequency of 2oo 2 ) excludes the influence of the uniform 
magnetic field generated by the current -carrying string 
1 which forms an offset signal at the input of the 
amplifier A when the string 1 is unperturbated from its 

30 straight line configuration. This increases the dynamic 
range of the magnetic gradiometer. It should be noted- 
that the uniform magnetic field produced by the 
alternating current of frequency co 2 through the string 1 
is also variable with the frequency co 2 and for this 

35 reason is cut off by demodulation at a frequency of 2co 2 . 

The string 1 is pumped with a sinusoidal current 
I(t) = isin(co 2 t) from a signal generator 3. The same 
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current is used to form a reference signal at twice the 
pumping current frequency 2co 2 by means of a multiplier 5 
and a high-pass filter 7. The frequency co 2 of the 
pumping current is twice the fundamental frequency of 
5 the string. The reference signal and the output from 
the low-noise amplifier A are input to a synchronous 
detector 9. The output voltage V out of the synchronous 
detector is directly proportional to a magnetic gradient 
B x2 around the string 1 and is the output of the device. 

10 The string 1 and the pick-up coils LI, L2 are 

situated inside a cryogenic nitrogen bath 11 at 77K. 
This increases greatly the signal-to-noise ratio by 
minimising environmental thermal noise influence and 
Johnson noise inside the pick-up coils themselves, as 

15 well as increasing the stability of the mechanical 

parameters of the string 1, in particular the stability 
of the amplitude of its forced resonant oscillations. 

The prototype shown in Figure 1 was tested in 
laboratory conditions. The string was stretched inside 

2 0 a ceramic housing and a second resonant mode of 4 84 Hz 

was obtained. Two coreless pick-up coils of 350 turns 
each were used, coupled to a commercially available 
Stanford Research Lock- In amplifier having an input 
noise of approximately 6 nV per root Hz. The output 
25 signal was detected at the doubled 484 Hz frequency (968 
Hz) and the bandwidth was chosen to be 0.3 Hz. 

A small magnet was moved beside the experimental 
set-up, and a clear magnetic gradient signal was 
observed. During the experiment the string was shaken 

3 0 and moved without the magnet present. The apparatus did 

not respond to this movement. However, a DC response to 
the position of the magnet was produced when the magnet 
was returned . 

In another experiment represented schematically in 
35 Figure 2, the magnetic gradiometer prototype 13 of - 
Figure 1 was placed onto a mobile frame and moved 
continuously along a 2 meter long track 15 over a 
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magnetic field source 17 (shown in the Figure 2 by a 
dark spot) . The magnetic field source 17 generated a 
field of differential forces proportional to B x x I z 
represented schematically in Figure 2 by arrows 19. The 
5 approximate speed of the prototype was one centimetre 
per second. The typical sensor response was recorded 
from the screen of an oscilloscope and is shown in 
Figure 2 as a plot of the output voltage V out against 
time t. 

10 Thus, the experiments have clearly shown a 

differential response of the string sensor typical to 
magnetic gradiometry. By using a commercial HT C 
superconducting quantum interference device (SQUID) 
instead of the low-noise amplifier A, the sensitivity to 
15 magnetic gradients would increase to a level so high 

that it would be easy to detect, for example, treasure- 
bearing remains of a submarine or of a ship from a 
reasonable distance . 

Another embodiment of the invention, according to 
2 0 which a very high sensitivity to magnetic gradients can 
be achieved without using SQUIDs, has been tested 
experimental ly . 

In this second embodiment, the string 1 is pumped 
with an additional carrier- frequency current and 
25 inductively coupled to a resonant bridge comprising two 
pick-up coils connected in a gradiometric configuration 
in parallel with a capacitance and tuned to the carrier 
frequency. This frequency is chosen to be far removed 
from the string's effective mechanical bandwidth which 
30 is normally limited to a few kHz even for very stiff 

strings. In this case, mechanical displacements of the 
string give a low-frequency modulation (at the rate of 
the second natural mode of the string) of the carrier- 
frequency voltage across the bridge. 
35 The string and the resonant bridge are cooled to 77 

K to reduce the effect of thermal noise on measurements. 
A double-lock-in scheme provides, firstly 
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amplification and detection of the carrier signal and, 
secondly, amplification and detection of the low- 
frequency envelope with an amplitude proportional to the 
measured magnetic gradient. 
5 The prototype measuring system according to the 

invention was tested in a gradient field and the results 
of the test are shown in Figure 3 . The upper curve 
shows the magnetic gradient signal in the time domain. 
The amplitude of the magnetic gradient signal is 

10 approximately 16 nT/m, the frequency of the signal is 
0.3 Hz, and the scale of the screen is 10 seconds per 
division. The lower curve shows the spectrum of the 
magnetic gradiometer output signal. The signal 
frequency is marked by the cursor. The scale is 0.2 Hz 

15 per horizontal division and 10.0 dBm per vertical 

division. The estimate of the signal to noise ratio 
gives a sensitivity threshold of approximately 0.36 nT/m 
per root Hz, for this particular sensor configuration. 
The observed noise is not instrumental but reflects the 

20 magnetic environment around the experimental set-up. 
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Claims 

1. Apparatus for the measurement of magnetic fields 
comprising : 

5 a string fixed at both ends, held under tension and 

arranged to carry a current ; 

and sensing means arranged to detect the transverse 
displacement of said string from an unperturbated 
position, due to a magnetic field acting on said string. 

0 

2. Apparatus as claimed in claim 1, wherein said 
sensing means comprises at least two sensors 
symmetrically longitudinally positioned about the mid- 
point of said string. 

5 

3. Apparatus as claimed in claim 1 or 2 , wherein the 
string is arranged to carry an alternating current. 

4. Apparatus as claimed in claim 3, wherein the 

0 frequency of the alternating current is substantially 

equal to the frequency of a natural ant i -symmetric mode 
of said string. 

5. Apparatus as claimed in claim 4, wherein the 
5 frequency of the alternating current is twice the 

fundamental frequency of the string. 

6. Apparatus as claimed in any of claims 3 to 5, 
wherein the sensing means is arranged to generate a 

0 detection signal due to the transverse displacement of 
the string and the apparatus further comprises 
demodulating means arranged to demodulate the detection 
signal with a demodulation signal, the demodulation 
signal having a frequency twice that of the alternating 

5 current . 
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7. Apparatus as claimed in any preceding claim, 
wherein said sensing means comprises at least one pick- 
up coil arranged such that, in use, a current is induced 
in the coil by the current passing through the string. 

5 

8 . Apparatus as claimed in any preceding claim, 
wherein said sensing means comprises means for detecting 
the transverse displacement of said string in two non- 
parallel planes. 

10 

9. A method of measuring magnetic fields comprising: 
providing a string fixed at both ends,, held under 

tension and carrying a current; and 

detecting the transverse displacement of said 
15 string from an unperturbated position, due to a magnetic 
field acting on said string. 

10. A method as claimed in claim 9, wherein the string 
carries an alternating current. 

20 

11. A method as claimed in claim 10, wherein the 
frequency of the alternating current is substantially 
equal to the frequency of a natural ant i -symmetric mode 
of said string. 

25 

12. A method as claimed in claim 11, wherein the 
frequency of the alternating current is twice the 
fundamental frequency of the string. 

30 13 . A method as claimed in any of claims 10 to 12, 
wherein the step of detecting the transverse 
displacement of the string includes demodulating a 
signal due to the displacement of the string with a 
demodulating frequency which is twice the frequency of 

35 the alternating current through the string. 
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14 . A method as claimed in 
wherein the displacement of 
electromagnet ically . 

15. A method as claimed in 
wherein the displacement of 
two non-parallel planes. 



PCT/EPOO/04270 

15 - 

any of claims 9 to 13, 
the string is detected 

any of claims 9 to 14, 
the string is detected in 
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